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Abstract 

In the 1 = sector there are more scalar mesons than can fit in one qq nonet. 
Consequently, many have claimed that there is in fact more than one multiplet, per- 
haps both qq and qqqq. Such proposals require the existence of at least two strange 
isodoublets (and their antiparticles) . The current PDG Tables list just one state, the 
-fCg(1430), while fits to data with Breit-Wigner forms and variable backgrounds can 
accommodate a k(900) too. Whether a state exists in the spectrum of hadrons is not a 
matter of ability to fit data along the real energy axis, but is completely specified by the 
number of poles in the complex energy plane. Here we perform as model-independent 
an analytic continuation of the LASS irK scattering results between 825 MeV and 
2 GeV as presently possible to determine the number and position of resonance poles. 
We find that there is a i^ 1430 )' but no k(900). The LASS data cannot rule out the 
possibility of a very low mass k well below 825 MeV. 

PACS numbers: 11.55.-m, 14.40. Ev 

Keywords: Analytic Properties, Scalar Mesons, Strange Mesons. 



1 Introduction 

For many years the scalar mesons have caused controversy. There is little consensus on the 
composition of many states, whilst for some their properties and even existence is a matter 
for debate. In the latest edition of the PDG Tables |IJ, there are four scalar-isoscalars below 
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1.55 GeV, with a fifth at 1.7 GeV. This is obviously too many for a standard qq nonet. 
However, many QCD-motivated models predict the existence of non-gg mesons, such as qqqq 
states 0, KK molecules [[| and glueballs [§]], and it is precisely in the isoscalar sector 
that these unconventional mesons are most likely to be found. This excess of isoscalars, 
together with the existence of two isovectors, the a (980) and a (1450), has been suggested 
as evidence for two scalar nonets 0, |J: the conventional qq nonet lying around 1.4 GeV 
with an unconventional one centred around 1 GeV. However, the PDG lists only one pair of 
scalar-isospinors in this mass region, the _K"q(1430). This has led some authors to postulate 
a light strange meson, known as the k p], [7], § . Evidence for this resonance has been claimed 
within certain models || [it], [U], [12| , whilst other studies dispute this |L3|, Fits to 

ttK production data, using Breit-Wigner forms, some with arbitrarily varied backgrounds, 
others with theoretically motivated backgrounds, allow a k to be accommodated. However, 
the existence of a state is not a matter merely of the quality of fit to data along the real 
energy axis. Indeed, a state is wholly specified by there being a pole of the S— matrix in the 
complex energy plane on the nearby unphysical sheet. Such specification, of course, does 
require experimental information in the relevant region of the real energy axis. 

Where this light strange meson has a mass around 900 MeV, it comes within the energy 
domain testable with the highest statistics results on ttK scattering, which come from the 
LASS experiment |15|. In this paper we present the results of a close to model independent 
determination of the number and position of poles in the I = ~, J = ttK scattering 
amplitude based on these LASS data above 825 MeV. The technique relies wholly on the 
analytic properties of S— matrix elements and requires no artificial separation into resonance 
and background contributions. We also consider continuations using results from an earlier 
SLAC experiment [TE\. In Sect. 2, we present the method and then in Sect. 3 describe the 
experimental inputs we use. In Sect. 4, we apply the analytic continuation methodology first 
to model data to assess its capabilities and limitations. In Sect. 5, we go on to study the 
analytic continuation using the real experimental results. We find there is no k(900). Since 
the method is one of continuing experimental data and the LASS experiment only provides 
information above 825 MeV we cannot rule out the possibility of an even lower mass k. In 
Sect. 6, we conclude. 
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2 Method 



The method we use is due to Nogova et al. [|17[]. It combines simple statistics with the 



analytic properties of the scattering amplitude to locate the positions of poles. As usual, we 
expect resonances to correspond to poles on unphysical sheets |18|] . 

2.1 Mapping 

In order to locate poles in the complex plane from scattering data on partial wave amplitudes 
along the real axis, we must perform an analytic continuation. To maximise the domain in 
which this continuation is valid, we begin by conformally mapping the partial wave amplitude 
into the unit disc. As shown in Fig. [I], for ttK scattering the s— plane is split into two regions 
and in this case we map the region outside the circular cut. The mapping is designed so that 
the cuts of the partial wave amplitude in the s— plane are mapped onto the circumference of 
the circle in the z— plane. This allows us to continue the amplitude throughout the complex 
z— plane. 
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Figure 1: The cut structure of the tcK partial wave amplitude, where 
Sth = (% + rn n ) 2 and the radius of the circular cut is s c = m 2 K — m^. 



The mapping is accomplished in two steps. Firstly, 




(1) 
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which maps the real cuts in the s— plane onto the positive real axis of the y— plane and the 
circular cut onto the negative real axis. Then to map the twice cut y— plane into the unit 
disc, we define 

z( S ) wyfya-y/vw-v^) , (2) 

ipvy( s ) + Vv( s ) -y( s th) 

where (3 is a real parameter which is chosen so that the region of interest in the s— plane 
is mapped close to the imaginary axis in the z— plane, thereby minimising the distance the 
continuation must cover. The points s = s t h and s = s c are fixed for any value of /3, being 
mapped to z = 1 and z = —1 respectively. Notice that the cuts in the upper half of the 
s— plane (such as the region where physical data lie) are mapped to the upper semi-circle in 
Fig. 0. 

It is obvious from this mapping that the physical region in the s— plane only covers a 
fraction of the circle. Even if we had data over an infinite range of energies, we could never 
complete the circle, see Fig. ^. Nevertheless, we only need analytically continue very close 
to the region of data (i.e. the solid arc in Fig. We shall return to this point later. 

From Fig. [| we can see that the mapping is highly non-linear in the following sense. As 
we increase the energy from threshold to infinity, the proportion of the circle covered by each 
increment falls very sharply. It is clear that the region between threshold (at 633 MeV) and 
the start of the data (at 825 MeV) covers a much longer arc than the region in which we 
have data (825 MeV to 2.51 GeV) and the region between the end of the data and infinite 
energy is much smaller still. This compression of the high energy region would suggest that 
the method is less sensitive to higher mass resonances. This apparent weakness actually has 
some benefit. As we explain in the next section, we test the scattering amplitude against the 
hypothesis that it contains a given number of poles. The computational burden increases 
rapidly as we increase the number of poles, so if we had to account for all the possible 
radial excitations in a given channel the calculation would soon become computationally 
prohibitive. Consequently, the procedure is only practicable for the lowest few states. 
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Figure 2: The z-plane showing how points in s-plane map. The mapping 
parameter (3 is chosen so that s = (1.4 + 0. 15i) 2 GeV 2 is mapped to the imag- 
inary axis. The thicker line shows the arc covered by the LASS data. The 
inset shows an enlargement of the key region close to 2 = i, from where we 
analytically continue. The symbols mark particular values of s as follows: 
♦ s = s th , □ s = oo, Vs = -s c , o S = is c , s = -j^(s c + is c ), • s = s c . 



2.2 Analytic Continuation 

For now we assume that we have a scattering amplitude, with errors, Aj,Q defined at discrete 

points all the way around our circle \z\ = 1. We define \8Zi\ as the average distance between 

the ith data point and its two nearest neighbours in the complex z— plane shown in Fig. £| 

In regions where these discrete points are densely packed, our scattering amplitude is most 

tightly controlled and so we weight the error on each point by the density of the data points 

in that region, e^iz) = A i y/\5zi\/27r. We make a smooth interpolation to the amplitude, 

Y(z), and the weighted errors, e(z), so that we have continuous functions defined on the 

entire circle. This weighting procedure ensures that the small region of the unit circle in 

Fig. |2], where we have experimental data, controls the analytic continuation into the nearby 

region where resonance poles are expected to sit. If F(z) is a square integrable function 

on the circle, C, then we can test how well this function describes the data through a x 2 > 

1 This method requires the errors on the real and imaginary parts of the amplitude to be equal. This 
simplification is arranged by taking the value of the error to be the larger of the error on the real part and 
the error on the imaginary part at any given energy. 
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denned by 



2 J_ 

X ' 2tt 



c 



F{z) - Y(z) 2 



dz\ . (3) 



We now introduce a non-zero weight function g(z), which is defined to be real analytic 
and constrained by \g(z) \ = e(z) around the circle. We expand the data and the trial function 
as Laurent series about the origin, so that 

y(z) = ^4 = E **** > = ~7T = £ G ^ • & 

Since partial wave amplitudes are real analytic, the coefficients and are real. Although 
we are expanding about the origin, the expansion is carried out round the circle, since this 
is where we have data. The singular Laurent coefficients will pick up any poles within the 
disc. Substituting Eq. (J|) into Eq. (|3|) gives 

oo oo 

X 2 = ( a ~ k ~ y - k "> 2 + ^( a k- Vkf ■ (5) 

k=l k=0 

The pole structure of our partial wave amplitude, Y(z), can then be determined by finding 
the test function F(z) which minimises the first summation in Eq. (|5|). 

If we want to test against the assumption that there are no poles in the data we must use 
a test function that is analytic, i.e. = for k > 0. Then, if the amplitude has no poles, 
the quantity Xo — 12k=i V-k wm ^ e zero 0- ^ we think that the scattering data has one 
pair of complex conjugate poles in the z— plane (as is the case where there is one resonance 
present) then we can write our test function as 

/(*) = + + Hz) 

Z — Zq Z — Z 

00 o?or„~fc-ii 



k=l 



E^+K») , (6) 



where h(z) is some analytic function. By comparing Eq. (|]) with Eq. we see that, in 

this case, the analytic continuation is carried out by setting a_fc = 2 KfazQ -1 ] an d then it 

is the quantity x\ — ^2k=i(V- k ~ 2 ^[a^Q -1 ]) 2 that has a x 2 distribution with iV degrees of 

freedom as long as \z \ <C 1. 

2 Actually, due to inevitable experimental noise, Xq wm n °t be exactly zero, but should fit a x 2 distribution 
with N degrees of freedom. 
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However, in practice \zq\ is close to one and the above procedure becomes unreliable. 
Consequently, we adopt the alternative approach of cancelling any pole in the data explicitly 
by the introduction of the so-called Blaschke pole-killing factor. This avoids the need to make 
the expansion explicit in Eq. and so significantly improves the convergence. For a pair 
of complex conjugate poles the Blaschke factor looks like 

(z-Zq)(z-zZ) 

(1 - zz )(l - zz*) ■ {t) 

We define the function 

and then the quantity J2k=i(y-k) 2 wm indeed obey a x 2 distribution with N degrees of 
freedom, if the actual amplitude contains one pair of complex conjugate poles. 



3 Experimental Input 

Our experimental information on the S— wave 7rK partial wave amplitude comes in the form 
of the magnitudes, a(s), and phases, <f>(s), measured, for instance, by the LASS experiment. 
The partial wave amplitude on sheet I is normalised thus 

where p(s) = 2q/y/s and q, the cm. 3-momentum 



1 [s - (m K + m^) 2 }^ - (m K - m n ) 2 } 

9= 2V ; • (10) 

ttK scattering has two possible isospin channels, 1 = 1/2 and / = 3/2. The resonances 
we consider are isodoublets, as there are no known resonances in the I = 3/2 channel. 
This means that there is no need to separate the isospin components to determine the pole 
positions, whereas this is essential in a Breit-Wigner-type fit. We therefore carry out our 
analysis for the total S— wave data. However, as a check we also perform this for the 1 = 1/2 
results alone. Of course, the 1 = 1/2 data must be extracted from the total S— wave data by 
subtracting an assumed form for the 1 = 3/2 amplitude, thereby introducing an additional 
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layer of uncertainty. Moreover, this separation is only possible over a limited energy range. 
Nevertheless, the I = 1/2 component does provide a valuable check. 

As is well-known, resonance poles do not appear on the physical sheet in the energy plane. 
Consequently, we must move to the relevant unphysical sheet. For purely elastic scattering 
this would be sheet II. In practice for txK scattering, the rjK channel opens very weakly (in 
agreement with SU(3)p expectations) and so any inelasticity can be safely neglected, until 
one reaches the rj'K threshold. This channel opens in the region of a possible i^o(1430) and 
complicates the sheet structure. Any resonance would then be on sheet III0. We must take 
this threshold into account when we change sheets. The simplest way to pick the correct 
sheet is to change the sign of the phase, i.e. 

This moves us onto sheet II below the rj'K threshold and sheet III above it. This method 
of changing sheets is also valid for the total S— wave amplitude. We then map the data, 
including errors, as described above. 

As noted previously it is not possible to cover the circle completely with physical region 
data, since we do not have experimental results either close to threshold or on the circular 
and left hand cuts (see Fig. |1|). The exact proportion of the circle covered by the direct 
channel 7rK results depends on the choice of the mapping parameter (3. How we treat the 
near threshold arc of Fig. ^| is a little different for our model tests of Sect. 4 and when we 
consider the actual experimental data in Sect. 5. Consequently, we detail their treatment 
in the appropriate section. As we do not have information about what happens on the 
unphysical cuts, we simply make a guess. So that the guess does not unduly affect the 
results, we de-weight these guessed points by giving them very large errors and by ensuring 
that they are widely spaced (see Fig. With the top semi-circle now spanned by data 
points, we make an interpolation to give us a continuous function. We complete the circle 
by reflecting the interpolated data in the upper semi-circle onto the lower half so as to obey 
the Schwarz Reflection Principle. 

A suitable form for the weight function g(z), which fulfills all the conditions described 

3 The Riemann sheets are labelled by (ri,^), where the Tj are the signs of the imaginary parts of the 
complex phase spaces pj in particular channels, j = 1 for nK, j = 2 for rj' K. 
Sheet I is (+, +), sheet II is (— , +) and sheet III is (— , — ) 
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above is 




g(z) = exp > c n z n (12) 



where the c n are found from a Fourier cosine fit to lne(z). We take N = 100. The singular 
coefficients of the Laurent expansion of the data about the origin are just 

V-k = ± /^f \dz\ ■ (13) 
2vr J c g{z) 



4 Results from Model Data 

As a first application, the method described in Sect. 3 is tested on a model amplitude describ- 
ing the situation where a light, broad resonance, Ki, and a heavier, narrower resonance, 
are present^. The resonances are constructed using Jost functions, with the lighter, broader 
resonance being treated as the background at the second resonance, i.e. 

f I (s) = f((s) + e 2 ^fi(s) (14) 



on the first sheet, where 



fi( s ) = ^Ml (is) 

jW p(k + cj +idj)(k - cj +idj) ' 

c / \ ( 2kdi \ 

5i(s) = arctan I ^— - j , (16) 



k = -Vs-irriK + m^) 2 . (17) 

The four parameters {cj,dj} are chosen, so that the amplitude has poles at 
s t = (0.9±0.25i) 2 (GeV) 2 and s 2 = (1.4 ± 0.15i) 2 (GeV) 2 . Data are created for energies 



equivalent to the full LASS data-set fLo| , with the error on both the magnitude and phase 
fixed at 5%. Below the energy of the LASS data, we create a few equally spaced points 
between threshold and and 825 MeV. This will be discussed in more detail later. The am- 
plitude in the unphysical region is set to a real constant (incidentally equal to the amplitude 
at threshold) with very large errors (±5). These model data are illustrated in Fig. |3|. 
4 This generalises a model analysis discussed by Nogova et al. to two resonances. 
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Figure 3: Real and imaginary parts of f 11 as a function of 9, the angular posi- 
tion around the circle in degrees, for the model data-set. The circles represent 
data points, the solid line shows the interpolation to the data and the dashed 
line shows the interpolated weighted errors. 



These model data allowed us to verify that the Blaschke factor method is indeed more 
stable to variations of non-physical parameters (such as the number of terms in our definitions 
of the x 2 an d the weight function) and hence more accurate. The pole positions and x 2 ' s 
found by fitting over the first 40 singular coefficients using this procedure are shown in 
Table [[]. Also tabulated for comparison are the results obtained when the errors on the 
points in the unphysical region are halved. 

From Table [I], it is clear that the method is readily capable of identifying the number 
and position of poles even when the associated resonances are broad and overlapping. The 
decrease in x 2 when going from no resonance to one resonance is large, but far less significant 
than the fall when we ask is there one or two resonances. From this Table, one would have 
no doubt that there are two resonances in the model data and their positions are well 
determined. 

Evaluating the errors on the pole positions that we obtain using this method is never- 
theless not straightforward. The experimental errors are folded into the actual calculation 
and so their propagation cannot easily be followed. Moreover, as pointed out by Nogova 
et al. JTTj, the x 2 ' s obtained are not strictly statistical, so the standard confidence level 
techniques are not appropriate. Thus we resort to an order of magnitude estimate of the 
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Option 


No. of 






x 2 




resonances 


(r 6) 


(MeV) 






n 






4081 


1 


1 


(0.988,89.31°) 


1201 ± 131i 


281 




2 


(0.993,89.98°) 


1396 ± U2i 


0.5 






(0.943,88.64°) 


903 ± 234i 













3759 


2 


1 


(0.989,89.13°) 


1173 ± lllz 


417 




2 


(0.993,89.97°) 


1392 ± 140i 


1.0 






(0.942,88.61°) 


900 ± 233i 





Table 1: Pole positions and x 2 's for model data. Option 1 has the unphysical errors 
set to 5. Option 2 has the unphysical errors set to 2.5. 

errors obtained by varying the inputs that introduce uncertainty into the calculation. These 
inputs include the number of terms in the Fourier expansion defining the weight function, the 
number of terms in the summation used to evaluate \ 2 an d the treatment of the unphysical 
region. The uncertainty in the pole position found in the z-plane due to these changes in 
parameters is shown in Fig. |4], from which we see that the typical uncertainty in the z-plane 
is of the order of (4 + 20i) x 10~ 5 . 



Max. 
Energy 


xl 


xl 


xl 


A 

(MeV) 


Ai A2 
(MeV) 


22.51 


2320 


43.3 


2.02 


1395 ± 53i 


1938 ± 234i 1453 ± 137z 


12.51 


2361 


37.0 


1.98 


1394 ± 56i 


1951 ± 253i 1456 ± 133z 


8.51 


2389 


35.0 


1.96 


1394 ± 58i 


1959 ± 264i 1457 ± 130i 


6.51 


2380 


18.9 


1.84 


1395 ± mi 


2024±331i 1458±114z 


3.51 


2207 


11.4 


1.77 


1396 ± 78i 


2099 ± Uli 1454 ± 102i 


2.91 


1959 


7.33 


1.70 


1398 ± 84i 


2178 ± 788i 1445 ± 91i 


2.51 


1529 


4.58 


1.57 


1399 ± 92i 


899 ± 1457i 1422 ± 82i 



Table 2: % 2 's and pole positions for trial data described in the text. Maximum 
energies are in GeV. A is the pole position found when searching for one resonance. 
Ai and A2 are the pole positions found when searching for two resonances. Xo> x\-> 
and xl are t ne x 2 's assuming no, one and two resonances respectively. The actual 
pole positions are (1421 ± 119i) MeV for the k 2 and (1957 ± 106i) MeV for the « 3 . 

Of course the physical amplitude will contain many resonances with masses greater than 
those of our K\ and k^. As was stated earlier, this method is expected to be less sensitive to 
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Figure 4: Expanded view of the z-plane showing the spread of pole positions 
found by varying the input parameters. ♦: standard values. 0: double the 
number of terms in g(z). +: halve the number of terms in g{z). *: Double 
the number of terms in \ 2 . □: Halve the number of terms in \ 2 . o: Double 
the value of f(z) in the unphysical region. A: Halve the value of f(z) in the 
unphysical region. ■: Halve the errors in the unphysical region. •: Double 
the number of points in the unphysical region, i.e not in the direct channel. 



these poles. To see to what extent this is true a further test was carried out using trial data, 
up to various energies, containing an effective range background, a k 2 and another state at 
1.95 GeV, which we refer to as the K3. Data points were created at the same energies as 
the LASS data and above that at 40 MeV intervals with an error of 5%. The x 2 's and pole 
positions found are shown in Table |2|. From these results we can clearly see that, as the 
maximum energy of the data falls, the K3 becomes less necessary to describe the data in 
the z— plane. With data up to 22.5 GeV, there is a strong case for claiming that there are 
both the k 2 and the k 3 , but their parameters are not accurately found. When the data only 
extends to 2.5 GeV, the ratio of \\ t° x\ is small enough that one would not be confident 
in claiming that the data exhibited more than one resonance. While the position of the k 2 
is more stable, the parameters found for the K3 bear no resemblance to the correct values. 
Meanwhile, the pole parameters found when testing for just one pole are always recognisable 
as the k 2 . If we introduce realistic experimental errors and inelasticity, then the K3 becomes 
even more hidden. 
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Due to the non-linearity of the mapping a fixed absolute uncertainty in the z-plane will 
give different uncertainties in the .E-plane depending on the actual position. The minimum 
uncertainty in the position of the k 2 is of the order of 3 MeV on the real part and 5 MeV 
on the imaginary. For the k 3 the minimum uncertainty would be 5 MeV on the real part 
and 15 MeV on the imaginary part. In contrast, for the K\ the minimum uncertainty is of 
the order 1 MeV on both parts. The pole positions listed in Table |1] would suggest that 
these minimum uncertainties underestimate the true uncertainty and we can expect to do 
no better than an accuracy of 10 MeV on the real part and 20 MeV on the imaginary part 
of the pole position. 



5 Results from Real Data 




82 84 86 88 90 92 82 84 86 88 90 92 



Figure 5: Real and imaginary parts of f n as a function of 9 for the LASS 
S— wave data-set. Only the actual data points are shown, with unweighted 
errors. 



Having assessed the capabilities and limitations of the method using model data, we now 



apply it to the results from real experiments. From LASS [T3], we have the magnitude and 
phase of the tt + K~ S— wave amplitude from 825 MeV up to 2.51 GeV as shown in Fig. [|. 
Due to a Barrelet ambiguity, the LASS group find two partial wave amplitude solutions 
which only differ above 1.84 GeV. For our calculation we use their Solution A, but this 
choice of solution does not affect our conclusions regarding the k(900) and the ^(1430). 
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The total S— wave amplitude for ir + K — > tt + K is related to the amplitudes with definite 
isospinQ by 

fs(s) = / (1/2) + \ / (3/2) • (18) 

Since resonances are only expected in the I = 1/2 channel, it is natural to consider 
the effect of separating out this component. Such a separation requires a modelling of the 
1 = 3/2 contribution. Below 1.58 GeV, the LASS group provide a model of this contribution 
based on the parametrisation of Estabrooks fT9f . This allows us to apply the method used 
in Sect. 4 to the 1 = 1/2 amplitude alone, but only below 1.58 GeV. As a control on this 
procedure, we also consider the full S— wave data over this reduced energy range, which we 
refer to as the "Short S— wave data-set" . 

The effective range type formula provided by LASS gives one possible extrapolation of 
their data to threshold. This serves as a guide (and only a guide) to possible data points 
between threshold and the start of the data. Note that the superscript (/) is an isospin 
label. All equations for amplitudes in this section refer to the physical sheet. Their formula 
is 

fi 1 ) = | - (19) 

p (cot 4g - *) p ( cot 6bw ~ 

where the resonance term only appears in the 1 = 1/2 case and 

gcoUj& = Jo+^ ( V > (2°) 

cot d B w = o~F • ( 21 ) 

q is given by Eq. fllOD , and q r is its value at s = m 2 ,. The LASS fit^gave a/ 1 / 2 ) = 2.19 GeV -1 , 

b (i/2) = 3 U Gev-i, a WV = _i.03 GeV" 1 , = -0.94 GeV -1 , m r = 1.412 GeV and 

T r = 0.294 GeV. It is important to recognise that though the LASS effective range fit 

assumes the tail of a Breit-Wigner-like pole, this does not prejudge that such a pole exists 

in our analysis. We call this treatment of the low energy region Case A. 

5 This being a physical process, the imaginary part of each partial wave amplitude must be positive. 
The authors of Ref. |^o| have noted (after the present work was complete) that this favours solution A of 
LASS ||. 

6 These parameters were provided by W. Dunwoodie and are not the values quoted in the original paper. 
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While the above Case A is largely experimentally motivated, a more theoretically well- 
founded guide to low energy meson-meson scattering is provided by Chiral Perturbation 
Theory (%PT). This makes predictions for near threshold nK scattering that we can input 
into our analysis. However, the effect of the higher order corrections becomes larger and less 
immediately predictable as one goes much above threshold. Consequently, we only generate 
data points based on xPT within 100 MeV of threshold, with a precision encompassing the 



range of present calculations [pi] , |22| , f20f . We call this Case B. We let the method determine 
the interpolation between the low energy results and where LASS data begin. This 
avoids the need for us to enter into any debate about whether higher order corrections in 



%PT are summed better by the Inverse Amplitude method |22|, |23| or by explicitly including 
Such differences are large at 825 MeV. It is in this sense that we describe 



resonances 
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our method as "model independent" . 

In implementing either Case A or Case B, we only create a few low energy "data" points 
so as not overly to prejudice the results. These two alternative sets of "data" are shown 
in Fig. |6] together with the LASS experimental results below 1 GeV. In both Cases the 
unphysical cuts are treated as for the model data-set of Sect. 4. 





0.65 0.7 0.75 0.8 0.85 0.9 0.95 

Energy (GeV) 



0.65 0.7 0.75 0.8 0.85 0.9 0.95 

Energy (GeV) 



Figure 6: Magnitude, a s , and phase, 4> s (in degrees), of the n + K~ S 1 — wave 
below 1 GeV from LASS data (•), contrasting the near threshold fit by the 
experimental group (□) with the predictions of ;\PT (A). 

In Table |3|, we present the results of our analysis for the data shown in Fig. |5] and the two 
extrapolations shown in Fig. |6|. The Table shows the pole positions found for the full S— wave 
data, with full and halved errors on the unphysical points and for the Case A and Case B 
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Case 




Nn nf 

resonances 


(r,V) 


/« 

Y °pole 

^ivie v j 


A/ 2 













1373 




1 


1 


(0.994, 90.07 ) 


1433 ± 149? 


4.7 






2 


/ r\ c\c\r' c\r\ -i rn\ 

(0.995, 90.15 ) 


"1 A OO 1 "1 A O * 

1432 ± 148« 


1.8 


A 

A 






/ r\ ao'7 oo nno\ 

(0.987, 82.02 ) 


805 ± 13« 















1629 




2 


1 


(0.993, 90.06 ) 


1423 ± 157z 


9.0 






2 


(0.993, 90.06 ) 
(0.998, 86.98 ) 


1 mo I "1 r* * 

1423 ± 154z 
969 ± 6z 


8.8 













1340 




1 


1 


(0.994, 90.10 ) 


-\ A A A 1 1 r - 1 ' 

1444 ± 151z 


16.9 






2 


(0.994,90.07°) 


1436 ± 133z 


1.0 


B 






(0.854,95.14°) 


609 ± 238i 















1694 




2 


1 


(0.993,90.08°) 


1428 ± 168* 


33.2 






2 


(0.994,90.05°) 
(0.870,93.35°) 


1429 ± 132i 
645 ± 260i 


1.0 



Table 3: Pole positions and x 2 ' s f° r LASS S— wave data [0| from 0.825 GeV to 
2.51 GeV. Cases A and B are explained in the text and shown in Fig. |6]. 
Option 1: The amplitude in the unphysical region equals the amplitude at threshold, 
with an error of 5. 

Option 2: The amplitude in the unphysical region equals the amplitude at threshold, 
with an error of 2.5. 



continuations to threshold. It is worth pointing out that our method will 'find' exactly 
as many poles as it is asked to. So if we search for two resonances then positions for two 
resonances will be given no matter how many resonances are present in the data. Resonances 
that really are present in the data will be stable to changes in unphysical parameters and 
result in sizeable falls in the y 2 . Conversely, if the x 2 does not fall significantly between one 
and two resonances then we can conclude that the second resonance 'found' is not really 
present and the pole position given is meaningless. Likewise, any resonance found whose 
pole position changes wildly with variations in parameters not determined explicitly by the 
experimental ttK data used will also be an artifact. 

In order to corroborate the results shown in Table |]the results for the 1 = 1/2 and the 
shorter S— wave data-sets, using the Case A continuation to threshold are shown in Tables |] 
and ||. These Tables show a remarkable consistency, as one would hope for effects that are 
real. 
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wp LIU 11 


No of 

1NU. Ul 


^pole 


y °pole 


V 2 
A 




I cbOIldXLCtib 


(r R~\ 


I IVie V 1 













210 


1 


1 


(11.992,90.22 J 


1467 ± 224i 


2.1 




2 


(0.992,90.22°) 


1464±221i 


1.1 






(0.993,82.14°) 


808 ± 7i 













367 


2 


1 


(0.991,90.23°) 


1455 ± 245i 


5.3 




2 


(0.991,90.19°) 


1439 ± 237i 


2.4 






(0.911,100.1°) 


503 ± 238i 





Table 4: Pole positions and x 2 's for LASS / = \ data JT5] . 

Option 1 has unphysical errors set to 5. Option 2 has unphysical errors set to 2.5. 



Option 


No. of 


Zpole 


\J Spole 


x 2 




resonances 


(r,6) 


(MeV) 













170 


1 


1 


(0.991,90.03°) 


1471 ± 242i 


2.2 




2 


(0.985,90.07°) 


1422 ± 553i 


1.9 






(0.997,90.08°) 


1457 ± 80i 













283 


2 


1 


(0.991,90.26°) 


1459 ± 264i 


5.6 




2 


(0.990,90.23°) 


1440 ± 263i 


2.4 






(0.886,101.8°) 


513 ± 208i 





Table 5: Pole positions and x 2 's for the short LASS S— wave data [|TjJ. 

Option 1 has unphysical errors set to 5. Option 2 has unphysical errors set to 2.5. 

As a further check of our results, the same technique was applied to the total S— wave 
data from an earlier experiment |jl~6 . These data extend closer to threshold, so it is only 
necessary to create a point at threshold. Barrelet ambiguities lead to four possible partial- 
wave solutions, which differ only above 1500 MeV, so the choice of solution should not affect 
our conclusions regarding the k. We have carried out our analysis using solutions A and B 
and found the results to be qualitatively similar and so we show only those for solution B. 
The data are shown below 1 GeV in Fig. [7], which should be compared with the low energy 
treatment of the LASS data-set (see Fig. H). This new Case A is defined using an effective 
range formula, see Eq. [20|, with a/ 1 / 2 ) = 2.39 GeV -1 and c^ 3 / 2 ) = —1.00 GeV -1 . No resonance 
was assumed in this fit, so any lingering doubts that assuming the tail of the -K"q(1430) in 
the low energy extrapolation introduces a prejudice is clearly not there in this case. We also 
consider Case B using the scattering length predictions of x?T, as before. The results for 
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0.65 0.7 0.75 0.8 0.85 0.9 0.95 
Energy (GeV) 




0.65 0.7 0.75 0.8 0.85 0.9 0.95 
Energy (GeV) 



Figure 7: Magnitudes, a s , and phases, <fi s (in degrees) S— wave nK scattering 



amplitude below 1 GeV from Estabrooks et al. [ 16 1, to be compared with Fig. 6 



Case 


Option 


No. of 

resonances 


Zpole 

(r,9) 


(MeV) 


x 2 













553 




1 


1 


(0.991,90.20°) 


1446 ± 232i 


16 






2 


(0.986,89.83°) 


1277±231i 


12 


A 






(1.017,90.46°) 


1344 ± 444i 















812 




2 


1 


(0.990,90.16°) 


1411 ± 250i 


14 






2 


(0.988,89.92°) 
(1.01,90.47°) 


1323 ± 221i 
1564 ± 316i 


14 













551 




1 


1 


(0.991,90.20°) 


1447 ± 231i 


17 






2 


(0.986,89.84°) 


1286 ± 227i 


14 


B 






(1.015,90.44°) 


1382 ± 423i 















816 




2 


1 


(0.990,90.16°) 


1411 ± 253i 


16 






2 


(0.991,89.81°) 
(1.001,89.89°) 


1324 ± 159i 
1395 ± 18z 


16 



Table 6: Pole positions and x 2 's for Estabrooks et al. [16] S— wave data. The Cases 
A and B are described in the text. 

Option 1 has unphysical errors set to 5. Option 2 has unphysical errors set to 2.5. 
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the Estabrooks et al. data-set are shown in Table || 

For the real ttK experimental results, Tables display a consistency in identifying 
a single resonance which looks very like the ^(1430). The fall in \ 2 m going from no 
resonance to one resonance is always sizeable, ranging from a factor of 35 to a factor of 
290. In contrast when going from one to two resonances, the \ 2 does not always fall by a 
significant amount. Recall importantly that the present analysis technique is always most 
sensitive to the lightest resonance and this is clearly the _K"q(1430). So whilst the _K"q(1430) 
is overwhelmingly apparent, a lighter, broad resonance in the LASS range above 825 MeV 
is not. When using the LASS data and Case B, where the low energy behaviour is given by 
xPT, there is a sizeable decrease in x 2 for two resonances, but the very light and broad one 
is below the experimental ranges studied by both LASS and Estabrooks at al. 

Halving the errors in the unphysical region should not greatly affect the pole positions 
obtained, and from Table [I] we can see that this is indeed the case. From Tables 0-|6| we see 
that for the assumption that there is only one resonance present, then halving the unphysical 
errors or changing the near threshold continuation, has a small effect. However, when we 
force the amplitude to have two resonances the effect is much more noticeable. Generally, 
the resonance which looks like -K"q(1430) tends to stay in a similar position, but the second 
one moves around wildly. This is what one would expect if this second resonance is merely 
an artifact of forcing the method to find two poles that are not really present in the energy 
range of the data that are being analytically continued into the complex plane. 

We mentioned earlier that we would expect the method to be less sensitive to high mass 
resonances. This is borne out by these results. The LASS group provide strong evidence for 
a state at 1.95 GeV. With data to 2.51 GeV we would expect to see this state showing up, 
but the compression of the high energy portion of the amplitude has diluted the effect of the 
Kq(1950) so much that it is not needed to describe the data in the complex z— plane. This is 
to be expected from our study using model data, described in Sect. 4. The mapping procedure 
severely limits our ability to find states above 1800 MeV in ttK scattering. However, below 
that, and particularly down towards the lower end of the experimental range, the method is 
totally reliable and quite unambiguous. Consequently, we have no doubt that 

i. there is a fQ(1430), 

ii. there is no k(900), 

in experimental 7rK scattering data. 
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6 Conclusions 



We have here applied a method of analytic continuation in as model independent way as is 
presently possible to the existing data on tt + K~ scattering between 825 MeV and 2.5 GeV. 
We find that there is only one scalar resonance between 825 and 1800 MeV, which is readily 
identified with the Xq(1430). Our procedure requires no assumptions about the amplitude 
being described by Breit-Wigner forms with any particular form of background. It directly 
counts the number of poles of the S— matrix on the nearby unphysical sheet. This method 
provides a rigorous test of whether the k(900) exists, and we conclude it does not. This is so 
whether we use the older Estabrooks et al. Jl6| data or that of the LASS collaboration |i~5 |. 
Clearly, the question of the existence of a k of mass and width that does not intrude into 
the energy range explored by experiment cannot be addressed. 

The fact that there is only one strange scalar between 800 and 1800 MeV has implications 
for quark models that postulate both a qq and a qqqq nonet in this region Jl, |24| . Our results 
of course do not depend on this modelling. However, we cannot refrain from commenting 
that we believe that other calculations clearly show the scalars are different from other states 
with underlying qq composition. Their decay modes, overwhelmingly to two mesons, mean 
that the physical hadrons do have a large qqqq component in their Fock space [[11], [1^, [13], ^5| . 
Thus for many purposes the scalars below 1800 MeV (particularly those with flavour) do 
behave like 4-quark states, even though they may well be seeded by single qq pairs at the 
bare level. However, such calculations and speculations are for elsewhere. 

We conclude that the data on S— wave irK scattering exhibit just one pair of complex 
conjugate poles between 0.8 and 1.8 GeV. This means that there is only one resonance 
present in this channel in the region accessed by the presently highest statistics experiment. 
This resonance has a mass around 1400 MeV and a width of about 300 MeV. The k(900) 
does not exist. 
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